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Abstract. We present circular polarization measurements of circumstellar H2O masers. The magnetic fields in 
circumstellar envelopes are generally examined by polarization observations of SiO and OH masers. SiO masers 
^\ ' probe the high temperature and density regime close to the central star. OH masers are found at much lower 

densities and temperatures, generally much further out in the circumstellar envelope. The circular polarization 
detected in the (616-623) rotational transition of the H2O maser can be attributed to Zeeman splitting in the 
^ ' intermediate temperature and density regime. The magnetic fields are derived using a general, LTE Zeeman 

analysis as well as a full radiative transfer method (non-LTE), which includes a treatment of all hyperfine com- 
ponents simultaneously as well as the effects of saturation and unequal populations of the magnetic substates. 
The differences and relevances of these interpretations are discussed extensively. We also address a non-Zeeman 
■ interpretation as the cause for the circular polarization, but this is found to be unlikely. We favor the non-LTE 

' analysis. The H2O masers are shown to be unsaturated, on the basis of their line widths and the lack of linear 

, polarization. The field strengths are compared with previous detections of the magnetic field on the SiO and OH 

masers. Assuming a r^^ dependence of the magnetic field on the distance to the star, similar to a solar-type 
magnetic field, our results seem to indicate that we are probing the highest density maser clumps at the inner 
edge of the H2O maser region. This allows us to estimate the density of the clumps, and the extent of the H2O 
Q ' maser region. We show that the magnetic pressure dominates the thermal pressure by a factor of 20 or more. We 

^ , also give an order of magnitude estimate of the magnetic field on the surface of the stars. In particular we discuss 

C/3 ■ the differences between Supergiants and Mira variable stars. 

a 
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1. Introduction the central star, at radii of 5-10 AU where the SiO maser 

emission occurs (e.g. Kemball & Diamond, 1997). The 

A large majority of stars go through a period of high mass ^g^sured circular polarization indicates magnetic field 

loss at the end of their evolution while climbing the asymp- gt^^j^gths of « 5-10 G, when assuming a standard Zeeman 

totk giant branch (AGB). This mass loss, of the order of interpretation. However, a non-Zeeman interpretation has 

10 to 10 Mo/yr, is the main source for replenishing ^^^^ proposed by Wiebe & Watson (1998), which only 

interstellar space with processed materials. Thus the mass ^^^^-^^^ ^^j^^ ^ 30 ^^^-^ ^^^^^ densities and 

loss mechanism is an important subject of study. In late temperatures and generally much further from the star, 
type stars, this high mass loss produces a circumstellar en- ^^^^^ observations measure fields of « 1 mG (e.g. 

velope (CSE) in which several maser species can be found, gzymczak & Cohen, 1997; Masheder et al., 1999) and there 
These masers, especially SiO, H^O and OH, are excellent ^-^^^^ ^-^^^^^ 2eeman origin of the polarization, 

tracers of the kinematics of the CSEs. The role of mag- ^^^.^ recently, no information on the magnetic fields at 

netic fields in the mass loss and formations of CSEs is still distances of a few hundred AU from the star was avail- 

unclear. Polarization observations of circumstellar masers ui t^i . . ,1 . 1 4-1, xj • 

able, ihis is the region where the H2O maser emission 
can reveal the strength and structure of the magnetic field o , ■ i- i 1 

... occurs. Because water is a non-paramagnetic molecule, 

throughout the CSE. Observations of SiO maser polariza- j , ■ j-u i- c u • • -c 4-1 

° ^ determination 01 the magnetic held is sigmfacantly more 

tion have shown highly ordered magnetic fields close to -i-rn u rpi r. , 1 i. ^ ■ -i-i, tj n, 

° ° dirfacult. i he faelds expected m the H2O maser region are 

~c J a ~1 T~r f i • 4- i -j • n Stronger than the fields measured for the OH mascrs; since 
bend ojjprint requests to: WV (vlemming©strw.leidenuniv.nlj ' 
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Fig. 1. Synthetic V-spectra for an external field of 50 mG 
and a line width of Aul = 0.8 km s~^. 
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Fig. 2. B — Py relations. The solid lines represent the 

F = 7 — 6 transition for different Awl. The short- and 
long-dashed lines are the F = 6 — 5 and the F = 5 — 4 
transition respectively. The dotted line is a typical best 
fitted model including all three hyperfine lines. 



they probably occur in gas that is a factor of 10 — 1000 
more dense and closer to the central star than the OH 
masers. So we expect fields of a few tens to a few hundred 
mG. For these fields the Zeeman splitting of H2O is ex- 
tremely small, only » 10~^ times the typical half-power 
width of the H2O maser line (A^l « 30 kHz). However, 
Vlemmings et al. 2001 (hereafter VOX) have shown that in 
the presence of such magnetic fields the Zeeman splitting 
can be detected with high spectral resolution polarization 



observations. They presented the first results of circular 
polarization measurements of circumstellar H2O masers 
on the masers around the supcrgiant star S Per. A 'ba- 
sic', LTE, Zeeman analysis was used to infer the magnetic 
field strength along the line of sight. This analysis did not 
include the interaction between the various hyperfine com- 
ponents of the 22 GHz H2O transition. The observations 
confirmed the expected field strength, finding a field of 
« 250 mG. The method used in VOl was similar to the 
method used by Fiebig & Glisten (1989, hereafter FG) to 
analyze the circular polarization of strong interstellar H2O 
maser features. Here we present observations for a sam- 
ple of 4 stars. We discuss the data reduction path, and 
in addition to the FG method, we fit our observation to 
the theoretical, non-LTE, models of Nedoluha & Watson 
(1992), hereafter NW. These include hyperfine interaction 
as well as saturation effect. The non-Zeeman interpreta- 
tions presented in Nedoluha & Watson (1990) and Wiebe 
& Watson (1998) are also discussed. Finally we present 
measurements of the linear polarization. 

2. Theoretical Framework 

2.1. Standard Zeeman interpretation 
2.1.1. Basic LTE analysis 

Our first treatment of the standard Zeeman interpreta- 
tion closely follows the analysis performed in FG and is 
more thoroughly discussed in both VOl and FG. This is a 
standard method that assumes LTE. Thus, the narrowing 
and rebroadening of the lines as a result of maser satura- 
tion are not reproduced. We have included the possibility 
of multiple masing hyperfine components, assuming the 
F = 7 — 6 hyperfine component to be dominant. Each hy- 
perfine component will split into 3 groups of lines (17+ , a~ 
and tt). For a magnetic field B parallel to the line of sight, 
the Zeeman pattern consists of the two circular polarized a 
components only. The right- and left-handed (RR and LL) 
spectra, corresponding to the cr^ components will only be 
slightly shifted against each other. This corresponds to a 
characteristic frequency shift of the order of Avz ~ 10^ 
Hz-B[Gauss]- As a result, the observed V-spectrum (RR- 
LL) will be a sine-shaped function, corresponding to the 
derivative /' of the total power spectrum. The amplitude 
of this function depends on the maser line width, the mag- 
netic field strength, and on which hyperfine components 
actually contribute to the maser. By calculating synthetic 
V-spectra wo find the following relation for the percentage 
of circular polarization: 

— (^^ax ^'^in)/-^max 

= 2 • Af_f' • B[Gauss]Cos6'/AvL[km s-1]. (1) 

Here Vmax and Vmin are the maximum and minimum of 
the synthetic V-spectrum fitted to the observations. Aul is 
the full width half-maximum (FWHM) of the total power 
spectrum, and /max is the peak flux of the maser feature. 
B is the magnetic field strength and 6 the angle between 
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Fig. 3. Fractional circular polarization in units of 0.1%. The solid line is the total intensity Stokes /, the long-dashed 
hue is Stokes V. Left has Wthormai = 0.6 km/s, right has 1.0 km/s. The models are labeled for emerging brightness 
temperature T^^Ail = 10^^, 10^° and 10^ from top to bottom. The solid vertical hne indicates the peak and the short 
dashed lines the FWHM of the profile. 



the magnetic field lines and the line of sight. The ^f-F' 
coefficient depends on the masing hyperfine components. 
The Ap-pr coefficients have been determined by calculat- 
ing Py from synthetic V-spectra, determined for a series of 
magnetic field strengths B, and for the different hyperfine 
components. Some examples of synthetic V-spectra for the 
three hyperfine lines are shown in Fig. |l|. The B - Py 
relation is shown in Fig. for different hyperfine compo- 



nents and line widths. For the different hyperfine transi- 
tions separately, FG find ^f-F' = 0.013, 0.08 and 0.01 for 
the 7 — 6,6 — 5 and 5 — 4 transitions respectively. For the 
fitted combinations of hyperfine components we generally 
find Af-f' ~ 0.011. 
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Fig. 4. Linear polarization for 9 = 75°. The dashed line is linear polarization is given in units of 10%. Labelling is 
similar as Fig. || 



2.1.2. Full radiative transfer, non-LTE analysis 

The transfer equations for the polarized radiation of the 
22 GHz H2O masers in the presence of a magnetic field 
were solved in NW. We have redone the calculations to 
fit the theoretical total intensity and circular polarization 
spectra to our observations. They were performed for the 
combination of the three strongest hyperfine components 
which have previously been verified to make the only sig- 
nificant contribution (Deguchi & Watson, 1986; Nedoluha 



& Watson, 1991). We solve a related set of 99 equations 
of state for all the magnetic substates of the relevant hy- 
perfine components. 

The goal of NW was to asses the validity of Eq. |^ 
when specific effects of maser radiation are included in 
the analysis. The LTE treatment described above for ex- 
ample, does not include: narrowing and rebroadening of 
the maser line with increasing optical depth, hyperfine 
interaction, unequal populations of the magnetic sub- 
states and the generation of circular polarization due to 
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Fig. 5. a) Spectral line width versus emerging brightness temperature TbAJ7. The solid, dotted and dashed lines have 
A?;thermai = 0.6, 1.0 and 2.0 km/s respectively, b) Py versus TbAfi. Labeling is similar to a), thin hues have 9 = 0°, 
thick lines correspond to 6* = 75°. 



changes in direction of the linear polarization. To solve 
the equations of state, some assumptions have been made. 
The Zeeman frequency shift g^l^ is assumed to be much 
greater than the rate for stimulated emission i?, decay 
rate T and cross- relaxation rate T^. In this regime the off- 
diagonal elements of the density matrix which describes 
the molecular states are negligible, which greatly simpli- 
fies the calculations. The maser is assumed to propagate 
nearly one-dimensionally. Additionally, the pumping has 
a Maxwellian velocity dependence and is the same for all 
magnetic substates. 

For expected magnetic fields of « 100 mG the val- 
ues of for the stronger hyperfine transitions are about 
1000 s~^. The rate for stimulated emission can be esti- 
mated by: 

R ~ AkThl^Q./'iTTnuo (2) 

Here Tb is the brightness temperature, AfJ the beam- 
ing solid angle for the maser radiation and A the Einstein 
coefficient. For TbArj = lO^^^ R w 100 s"^ Thus 
gVL 3> i?, required for ignoring the off-diagonal elements 
of the density matrix, is satisfied for most of the relevant 
brightness regime, up to TbAil = 10^^. For higher bright- 
ness temperatures the results for 9 ^ Q° may be unreliable. 
For = 0° the off-diagonal elements are not present, thus 
the results are reliable up to the highest values of TbAil. 

Detailed calculations of the cross-relaxation rate (T^), 
which is expected to be considerably larger than the de- 
cay rate (F), give Fi, = 2 s~^ at temperatures of 400 K 
(Anderson & Watson, 1993). So gfi > F^ > F is satis- 
fied as well. The calculations here use (F -I- Fi,) = 1 s~^, 
but NW have shown that the resulting intensities scale 



with the decay rate and the cross-relaxation rate as 
[fiux/(F + F,)]. 

To perform the numerical calculations, we had to chose 
the strength of the radiation incident onto the masing re- 
gion. The calculations presented here are based on a ra- 
diation of TbAi7 — 0.1 K sr. This continuum radiation is 
taken to be unpolarized. It was verified in NW that the 
results are insensitive to the chosen value, and we have 
confirmed this result. 

In contrast to the LTE analysis presented above, the 
results from the non-LTE transfer treatment are V-spectra 
which arc not anti-symmetric. This can be seen compar- 
ing the non-LTE spectra in Fig. || to the synthetic LTE 
spectra in Fig. |^. The spectra are not proportional to the 
derivative of the total power intensity profile, F, as used 
before, due to the blending of hyperfine lines with different 

gn. 

Fig. m also shows that for high TbAll, the total in- 
tensity spectra are increasingly less Gaussian in shape. 
The line shape thus enables us to estimate the saturation 
level of the maser. The model with TbAJ7 = 10^ corre- 
sponds to a completely unsaturated maser. The models 
with TbAll — 10^" and 10^^ correspond to slightly and 
fully saturated masers respectively. The non-LTE analy- 
sis also produces linear polarization for 9 ^ Q° . Stokes Q 
is shown in Fig. |[ which indicates that for 9 = 75°, linear 
polarization of up to k, 10% is observed for TbAil « 10^^. 
NW show that for slightly saturated masers, linear polar- 
ization of a few percent should be observed when 9 > 15°. 

The calculations have been performed for several dif- 
ferent values of intrinsic thermal width (wthormai)- This 
is the FWHM of the Maxwellian distribution of parti- 
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Fig. 6. Ap^F' coefficients as a function of TbAri. The 
lower panel shows 9 = 0°, the upper panel 6 = 75°. The 
solid, dotted and dashed lines correspond to Awthcrmai ~ 
0.6, 1.0 and 2.0 km/s respectively. 



cle velocities. Assuming a kinetic temperature of T K, 

^thermal « 0.5(T/100)l/2 (NW). 

The variation of line width with brightness tempera- 
ture is shown in Fig. |5^ for different values of f thermal. 
This relation does not depend on the angle 9. We see 
a gradual narrowing of the maser line, and rebroaden- 
ing occurs when the maser starts to become saturated, 
at TbAfi « 10^°. The relation between the fractional cir- 
cular polarization Py and Tb does depend on 9 and is 
shown for different values of 9 and tithormai in Fig-^- The 
fractional circular polarization decreases with increasing 
maser intensity as the total intensity grows faster than 
the circular polarization. Both these relation scale with 
[flux/(r -I- Fiy)], similar to the intensity proffies. 

From these relations we can now determine the coeffi- 
cients Ap-pi from Eq.^ which are shown in Fig.^. As can 
be seen, Af-f' depends strongly on wthcrmai. For small 
angles 9 there is no strong dependence on Tb, but for 
9 > 45°, Af-f' can change by over a factor of 6 for high 
Tb. Neufeld & Melnick (1990) argue that the low values 
of ^thermal Can be eliminated, on the basis of realistic H2O 
maser pumping schemes. Values as high as ^thermal ~ 2.0 
km/s indicate temperatures of more than 1000 K, which is 
expected to be too high according to the same analysis by 
Neufeld & Melnick. From linewidth analysis we conclude 
that ^thermal ~ 1-0 km/s is the most realistic estimate, 
which yields Af-f' ~ 0.018. This value does not decrease 
by more than 10 % if we take uthcrmai ~ 1.5 corresponding 
to the upper temperature limit of 1000 K. 



Fig. 7. 0-dependence of Eq. |lj for increasing emerging 
brightness temperature TbAil. The lines for TbAil = 10^ 
and 10^ coincide and are the same as the lines for lower 
brightness temperatures. 

The effect of the angle 9 has recently been examined 
by Watson & Wyld (2001) for the general case of a tran- 
sition with low angular momentum (J=l-0 and J=2-l). 
They conclude that for masers that are at least partially 
saturated the circular polarization does not decrease with 
increasing 9 until the magnetic field and the propagation 
direction are nearly perpendicular. We have performed 
similar calculations for the specific 22 GHz H2O maser 
transitions and the results are shown in Fig. |7| which con- 
firm for transitions with higher angular momentum, the 
analysis of Watson & Wyld. 

2.2. Non-Zeeman interpretation 

Circular polarization can also occur due to the change of 
the axis of symmetry for the molecular quantum states 
from parallel to the magnetic field, to parallel to the di- 
rection of propagation. When gft ^ R, the direction of the 
magnetic field is the axis of symmetry. When gH. ^ R, the 
molecules interact more strongly with the radiation than 
with the magnetic field, and the direction of propagation 
is the axis of symmetry. Initially gfl 3> R, but as the 
maser propagates through the medium and becomes more 
intense, gfl <C R could possibly be satisfied. This was dis- 
cussed by Nedoluha & Watson (1990) and was shown to 
lead to an overestimation of the magnetic field by up to 
a factor of 1000. However, the study in NW, and the line 
width observed here, show that the effect described above 
is unimportant for H2O masers. The line widths show that 
in almost all cases gfi 2> i? is satisfied. 
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Wiebe & Watson (1998) have shown that the propa- 
gation of Unear polarization can also result in circular po- 
larization. For linear polarizations of ~ 10%, the resulting 
circular polarization is of the same order of magnitude as 
the polarization due to the regular Zeeman interpretation, 
while magnetic fields could again be a factor of 1000 less. 
Below we indicate that also this non-Zeeman interpreta- 
tion is unlikely due to the lack of linear polarization. 

2.3. Fitting Method 

2.3.1. LTE Zeeman 

In the case of the LTE Zeeman models we first fit a com- 
bination of the 3 strongest hyperfine components to the 
total power spectrum. From this we get the line width Awl 
of the maser feature and the contribution to the line of the 
hyperfine components. For this combination we calculate 
a synthetic V-spectrum, such as shown in Fig. |l] for the 
separate hyperfine components. The synthetic spectrum 
is then fitted to the circular polarization spectrum. As 
our observations in VOl have shown that the observed V- 
spectrum can be narrower than the synthetic V-spectrum, 
we also allow for the narrowing of our synthetic spectra in 
the fit. Furthermore we need to remove the scaled down 
replicas of the total intensity. We fit the following function: 

Vobsi^v) = ai • V*{a2 ■ Av) + ■ /(Aw) (3) 

Here Av is the velocity in [km/s] measured from the 
peak of the total intensity, V* is the synthetic spectrum 
K;ynth- Wc fit for the parameters ai, 02 and a^. The param- 
eter ai controls the amplitude of the V-spectrum, while 02 
determines the amount of narrowing of the observed spec- 
tra with respect to the model spectra. Finally, the removal 
of the scaled down total intensity profile is determined by 
03. The amplitude ai of the best fitted model determines 
Knin and Vmax- Together with the peak intensity of the 
maser feature this gives Py. With the previously deter- 
mined line width Aul we calculate iJ| | using Eq. |^. 

2.3.2. Non-LTE 

For the non-LTE case, we fit out models to both total in- 
tensity and circular polarization simultaneously. We use 
Eq. ||with V* = T^nodci- Because the non-LTE models are 
already intrinsically more narrow, we first fix the width 
of the V-spectrum by setting 02 = 1.0 and only fitting ai 
and as. We have chosen to fit models with i^thermai = 0.8 
and 1.0 km/s, as the line widths of the maser features in- 
dicate that much higher intrinsic thermal widths are un- 
likely. We find that in some cases the circular polarization 
spectrum is still narrower than the model spectrum. For 
these features we allow for narrowing by releasing 02 . The 
combined fitting determines the circular polarization per- 
centage Py and the line width Awl. The coefficient ^f-f' 
is specifically determined for the best fitted model. Using 
this Eq. IT] again gives B\\. 



3. Observations 

The observations were performed at the NRAO]^ Very 
Long baseline Array (VLBA) on December 13th 1998. 
At the frequency of the 616 — 623 rotational transition of 
H2O, 22.235 GHz, the average beam width is ~ 0.5 x 0.5 
mas. This allows us to resolve the different H2O maser 
features in the CSE. The data were correlated twice, once 
with modest (7.8 kHz =0.1 km s~^) spectral resolution, 
which enabled us to generate all 4 polarization combina- 
tions (RR, LL, RL and LR). The second correlator run was 
performed with high spectral resolution (1.95 kHz = 0.027 
km s^^), necessary to detect the circular polarization sig- 
nature of the H2O Zeeman splitting, and therefore only 
contained the two polarization combinations RR and LL. 
We have performed 6 hours of observations per source- 
calibrator pair. The calibrator was observed for 1 1/2 
hours in a number of scans equally distributed over the 
6 hours. We used 2 filters (IFs) of 1 MHz width, which 
were overlapped to get a velocity coverage of « 22 km/s. 
This covers most of the velocity range of the H2O maser. 

3.1. Calibration of tiie Data 

The data analysis followed the method of Kemball, 
Diamond & Cotton (1995). The reduction path, shown in 
Fig. H, was performed in the AIPS data reduction pack- 
age. The first standard calibration steps were performed 
on the data-set with modest spectral resolution. We used 
the system temperature measurements provided with the 
data to perform the amplitude calibration for both the 
calibrators and the sources. Parallactic angle correction, 
flagging and single- and multi-band delay calibration were 
all done regularly on the calibrators observed with each 
source. During 1/3 of the observation time, the first IF 
suffered strong interference, which forced us to flag sev- 
eral of the frequency channels (« 0.1 MHz) in addition 
to the normal flagging. Also, for most of the observations 
IF 1 of the Los Alamos (LA) antenna was unusable. The 
solutions obtained at these calibration steps were copied 
and applied to the high spectral resolution data set. The 
complex bandpasses were then determined for both data- 
sets separately. Additional calibration steps were needed 
for accurate processing of polarization data. The gain ra- 
tio between the R- (right-) and L- (left-) hand polariza- 
tions was determined using the auto-correlation data of 
the reference antenna on a short scan of the maser source. 
This step contains the critical assumption of equal flux 
for the R- and L- handed circular polarizations. As a re- 
sult the circular polarization spectra are forced to be anti- 
symmetric, as discussed below. Later, during the analysis 
of the circular polarization spectra, we still have to cor- 
rect for small replicas of the total power profiles due to 
remaining small gain differences. After the gain calibra- 
tion, we used the calibrators to determine the R-L delay 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 
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Table 1. Star Sample 



Name 


Typo 


RA (J2000) 


Dec (J2000) 


Distance 


Period 


Vrad 


Peak flux 






m s ^ 




(pc) 


(days) 


(km/s) 


(Jy) 


S Per 


Supergiant 


02 22 51.72 


+58 35 11.4 


1610" 


822 


-38.1 


76.1 


U Her 


Mira 


16 25 47.4713 


+18 53 32.867 


189'' 


406 


-14.5 


12.5 


NML Cyg 


Supergiant 


20 46 25.7 


+40 06 56 


1220^= 


940 


0.0 


48.2 


VY CMa 


Supergiant 


07 22 58.3315 


-25 46 03.174 


1500"^ 


2000 


22.0 


244.1 



" Hipparcos,''van Langevelde et al. (2000), '^Dance et al.(26oi) , "Lada & Reid (1978) 



RR, LL, RL and LR 1 
7.8 kHz resolution | 








Sort and Flag data 
Parallactic Angle correction 



RR and LL only 
1.95 kHz resolution 



copy solutions 

— ► 



Velocity shift to LSR 
Reference Frame 
I 



Velocity shift to LSR 
Reference Frame 



PreHminary Bandpass conrection 
Amplitude calibration 
Single Band delay determination 
on calibrator 



Multi Band delay correction 
on calibrator only 



copy solutions 

— ► 



Complex Bandpass correction 



CompleK Bandpass correction 



R-L gain ratio determination 
on Maser source 



copy solutions 

— ► 



R-L delay offset and 
RR-LL phase offset determination 
on calibrator source 



copy solutions 

-► 



Fringe fit and Seif-Calibration 
on strong Maser feature 



Fringe Stand Sedf- Calibration 
on strong Maser feature 



Feed calibration on cross-power 
spectrum of the Maser source 



copy solutions 

-► 



Make I,Q,U image cubes 



Make I,V image cubes 



Fig. 8. The schematic data reduction path for om' spectral 
line polarization data 



offset and the RR-LL phase offset. The delay offset is ex- 
pected to be stable over the duration of the observation 
(Brown et al. 1989) and can be determined from a sin- 
gle continuum scan in which the cross-hand fringes RL 
and LR are well detected. The RR-LL phase offset varies 
over the observation and are determined from calibrators 
which are assumed to show no circular polarization. The 
solutions were again copied from the modest resolution 
date and applied to the high resolution data. Then fringe 
fitting and self-calibration for the two separate data-sets 
were performed on a strong maser feature. Finally, cor- 
rections were made for the instrumental feed polarization 



using a range of frequency channels on the maser source, 
in which the expected linear polarization or the frequency 
averaged sum of the linear polarization is low. After the 
solutions were applied to both data sets image cubes could 
be created. 

3.2. Sample 

We observed 4 late type stars, the supergiants S Per, 
VY CMa and NML Cyg and the Mira variable star U Her. 
They are listed in Table. |l] with type, position, distance, 
period and velocity. The sources were selected on the ba- 
sis of 2 criteria; strong H2O masers and previous SiO and 
OH maser polarization observations. The peak fluxes mea- 
sured in these observations are also shown in Table. |^. On 
the total intensity channel maps with high spectral resolu- 
tion, the noise is dominated by dynamic range effects and 
lies between « 0.08 — 0.3 Jy. On the circular polarization 
maps we have noise of « 0.01 — 0.03 Jy. 

The circular polarization of the SiO masers around 
VY CMa has previously been observed with a single dish 
by Barvainis et al. (1987). They find circular polariza- 
tion of 6.5 %, indicating a field strength of Bsio ~ 65 G. 
Observations of the 1612 MHz main line OH maser by 
Cohen et al. (1987) give Bqh « 2 mC. Smith et al.(2001) 
concluded from ground based and Hubble Space Telescope 
optical observation, that high magnetic fields of at least 
1 G are necessary to explain the outflow of matter ob- 
served in VY CMa. Cohen et al. also detected fields of 
« 2 mG on the 1612 MHz OH masers around NML Cyg, 
while Masheder et al. (1999) estimated a field of « 1 mG 
on the OH masers around S Per. On U Her, Palen & Fix 
(2000) observed pairs of R- and L-polarized maser fea- 
tures, for the 1665 and 1667 MHz OH masers. Although 
not many of these Zeeman pairs were found they estimate 
a magnetic field of Bqh ~ 1 mG. 

4. Results 

We have examined the strongest H2O maser features 
around the 4 stars observed. The results are shown in 
Table || for features with fluxes up to 10% of the bright- 
est maser spots. Some selected spectra are shown with 
fits of the LTE Zeeman models and the non-LTE models. 
When using the non-LTE models, fits with f thermal = 1.0 
km/s generally are best. As expected the splitting is small 
and the circular polarization is generally not more than 
0.5 — 1.0%. Magnetic field strengths obtained using the 
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Fig. 9. (left) Total intensity image of the H2O maser features around S Per (right). Total power (I) and V-spectra 
for selected maser features of S Per. The bottom panel shows the best fitting synthetic V-spectrum produced by the 
standard Zeeman interpretation (dashed line). The middle panel shows the best non-LTE model fit (dotted line). 
The corresponding total power fits are shown in the top panel. The V-spectra in the lower two panels are adjusted 
by removing a scaled down version of the total power spectrum determined in the fit by the 03 parameter in Eq.^ 
which is different for the LTE and non-LTE fits. The solid vertical lines show the expected position of minimum and 
maximum of the sine-shaped function in the general Zeeman case. 



LTE Zeeman models are given in column 7 while the re- 
sults obtained using the non-LTE models are shown in 
column 8. The features that show less than a Scr circu- 
lar polarization signal are considered non-detections. For 
these we have determined (absolute) upper limits using 
the ^F-F' coefficients obtained using the NW radiative 
transfer models. When using the LTE Zeeman analysis 
the limits increase by a factor of « L4. 

Features labeled * have V-spectra that did not allow 
accurate fitting, possibly due to blending effects. They do 
show circular polarization above the 3a level, although 
generally only slightly more. For these we have used Knin 
and Knax to estimate the magnetic field strength, and we 
assumed the best estimate for the coefficient Af-f' ~ 
0.018. The errors are the formal la errors. 

We have not been able to detect linear polarization in 
any of the 4 sources. This gives upper limits to the frac- 
tional linear polarization of « 1% for the weakest features 
down to « 0.01% for the strongest. 

4.1. S Per 

Fig. H shows the total intensity map of the H2O maser 
features around S Per, in which we can identify most of 
the features detected in earlier observations (Diamond ct 
al. 1987; Marvel 1997). Positions are shown relative to 
the feature labeled S Per b, discussed in more detail in 



vol. Circular polarization between 0.2 — 0.7% has been 
detected in 3 of the brightest features after a careful reex- 
amination of the data. This has also resulted in a slightly 
lower magnetic field strength value for S Per b than the 
one presented in VOl after optimizing the fitting routines. 
The figure also shows the total power and circular polar- 
ization spectra for 2 of the features, with fits for both the 
LTE and non-LTE Zeeman models. The velocity labeling 
has changed from VOl to correspond to the convention 
used by NW. We find the magnetic field pointing away 
from us on all features. From these results we estimate 
the magnetic field in the H2O maser region around S Per 
to be « 200 mG using the LTE Zeeman analysis or « 150 
mG using the non-LTE models. 

The maximum radius of the H2O maser region was 
determined to be « 60 mas by Diamond et al.(1987), cor- 
responding to « 100 AU. They conclude for a thick shell 
model that the inner radius of the water maser shell is 
« 30 mas, while the outer radius is w 70 mas. For com- 
parison, the OH maser shell extent was determined by 
them to be 80 x 30 mas, using MERLIN observations. 

4.2. VYCMa 

The total intensity map of the H2O masers around 
VY CMa is shown in Fig. ^ Circular polarization has 
been detected on 7 out of 10 of the brightest maser fea- 
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Table 2. Results 



Name 


Feature 


Flux (I) (Jy) 


V;.ad(km/s) 


AiiL (km/s) 


Pv (-10-") 


B cose (mG) LTE 


B cose (mG) non-LTE 


S Per 


a 


76.1 


-27.2 


0.88 






< 26 




b 


57.7 


-30.9 


0.77 


6.20 


-207±30 


-132±18 




c 


11.4 


-33.9 


1.21 






< 254 




d 


34.7 


-37.1 


0.54 


5.62 


-140±32 


-77±17 




e 


29.7 


-28.9 


0.61 






< 41 




f 


35.6 


-40.1 


0.53 






< 33 




g 


15.4 


-30.4 


0.92 






< 125 




h 


14.9 


-34.2 


0.48 






< 75 




i 


12.9 


-39.3 


0.69 






< 121 




• * 
J 


26.6 


-39.2 


0.64 


2.41 


-156±47 


-68±35 




k 


14.2 


-27.2 


0.51 






< 72 




1 


10.0 


-23.2 


0.44 






< 110 


VY CMa 


a* 


244.1 


10.8 


0.66 


3.70 


-116±21 


-151±21 




b 


155.5 


10.6 


0.67 


5.64 


182±30 


99±15 




c 


115.2 


12.9 


0.74 


2.73 


85±18 


112±19 




d 


47.9 


13.7 


0.83 


3.24 


129±31 


81±16 




e 


17.8 


15.0 


0.94 


4.10 


164±66 


166±55 




f 


58.5 


18.3 


0.80 






< 25 




g 


25.3 


27.1 


0.88 






< 55 




h 


43.1 


28.8 


0.61 






< 25 




i 


65.6 


25.2 


0.79 


2.23 


83±25 


48±14 




J 


126.2 


24.6 


0.60 


4.66 


133±26 


85±12 


NML Cyg 


a* 


48.2 


-21.2 


0.66 


5.9 


215±37 


118±15 




b 


10.5 


-20.4 


0.40 


13.8 


234±63 


221±41 




c 


5.42 


-20.2 


0.54 


19.2 


492±133 


336±67 




d 


3.66 


-20.0 


0.54 


21.6 


553±220 


428±149 


U Her 


a 


12.20 


-17.8 


0.75 


27.5 


858±161 


605±83 




b 


2.91 


-17.6 


0.50 


126.1 


2434±505 


1487±231 




c 


1.16 


-17.6 


0.73 






< 803 




d* 


2.46 


-19.2 


0.47 


49.6 


1181±297 


552±114 




e* 


2.07 


-19.2 


0.51 


41.1 


980 ±298 


726±189 




f 


1.36 


-19.3 


0.51 






< 337 



* See text 




MilliARC SEC 

Center at RA 07 22 58.29545 DEC -25 46 02.371 7 Av rkm/sl Av [km/sl 

Peak flux = 2.8461 E+02 JY/BEAM 

Leva = 2.846E+01 * (0.250, 0.500, 1, 3, 6, 9) 



Fig. 10. Similar to Fig. | for VY CMa 
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4.4. U Her 

The total intensity map for the H2O masers around U Her 
is shown in Fig. 13, U Her is the only Mira variable star 
in our sample and shows significantly weaker H2O masers. 
However as seen in the figure, we do detect circular po- 
larization of up to « 10%. The field strength in the H2O 
maser region of U Her is higher than those observed for 
the supergiant stars in our sample. We find strengths of 
« 2.5 G using the LTE Zeeman analysis or « 1.5 G using 
the non-LTE models. While the features U Her a and b 
show a clear sine-shaped V-spectrum, the features d and 
e have circular polarization signals only slightly outside 
the 3cr detection limit. This made fitting actual profiles 
impossible for those features. 

The size of the H2O maser region is found to be w 350 x 
200 mas (Bowers et al. 1989), which at 189 pc distance 
corresponds to « 65 x 40 AU. Chapman et al. (1994) find 
for the OH maser extent of U Her of 800 x 900 mas, ~ 
150 X 170 AU. 



Fig. 12. Total power (I) and circular polarization (V) of 
NML Cyg a. Dotted and dashed lines are the model of two 
slightly shifted spatially blended maser features. 

tures down to Py « 0.2%. Similar to the values found for 
S Per we find a magnetic field of « 200 mG using the LTE 
Zeeman analysis or « 175 mG using the NW models. The 
magnetic field points toward us for most of the sources. 
VY CMa a suffers from blending with nearby maser spots 
which made accurate fitting impossible, this can also ex- 
plain the change of sign. 

Diamond et al.(1987) give the radius of the H2O maser 
region for VY CMa to be « 150 mas, corresponding to 
PS 220 AU. The 1612 MHz OH maser radius is shown to 
be « 2400 AU in observations by Reid et al.(1981). 

4.3. NML Cyg 

We have only been able to detect a few H2O maser fea- 
tures around NML Cyg, which are shown in Fig. |l^. We 
find circular polarization from 0.6 — 2%. The spectra for 
NML Cyg b and c are shown in the same Figure. The 
V-spectrum for NML Cyg a is shown in Fig. |l^. The cir- 
cular polarization spectrum seems to indicate that this 
features actually consists of 2 heavily blended features of 
w 30 Jy. The magnetic field strength on this feature has 
been estimated using this simple model. Both the inter- 
pretations indicate the magnetic field strength to be « 500 
mG, higher than for both S Per and VY CMa. 

The H2O maser region extent is « 195 mas, corre- 
sponding to « 240 AU (Johnston et al. 1985). Diamond 
et al.(1984) give for the 1612 MHz OH masers a maximum 
extent of Ri 3 — 5 arcsec. 



5. Discussion 

5.1. Interpretation of the Polarization Results 

The observations show magnetic fields in the range from 
« 150 mG to « 1.5 G in the H2O maser regions of our 
sample of stars. The exact field strength depends on the 
interpretation used. Here we have used the LTE Zeeman 
analysis presented in EG and VOl, as well as a non-LTE 
model presented by NW. Generally, the results using the 
NW method give fields approximately 2/3 of the LTE 
Zeeman method. As shown above, the non-LTE models do 
not typically produce anti-symmetric shapes for the cir- 
cular polarization spectra. However, the anti-symmetric 
shapes observed are a necessary result of the data treat- 
ment, when it is assumed that the R- and L-polarization 
profiles are similar and when the scaled down replica of the 
total power spectrum is removed. We can thus fit both in- 
terpretations to the data by leaving the amplitude of the 
total power replicas as a free parameter as described in 
Section |2.3| . As a consequence, we are unable to directly 
distinguish between the different interpretations on the 
basis of the anti-symmetry of the V-spectra. 

The shape and line widths of the total intensity spectra 
of the different maser features indicate that for all stars 
the masers are unsaturated or only slightly saturated. As 
a result, even for the non-LTE case, we do not expect 
significant deviations from anti-symmetry. The line widths 
also indicate that ^thermal ~ 1.0 km/s is indeed the best 
estimate for the thermal line width. In the cases where we 
could not perform good fits to the spectra, we have thus 
used Af_f' — 0.018 as the best estimate. Since the masers 
are thought to be unsaturated the circular polarization 
decreases linearly with cos 9, as shown in Fig. 0. 

An important motivation for using the non-LTE ap- 
proach, is the fact that the circular polarization spectra 
show narrowing that is not expected in the LTE 
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Fig. 11. Similar to Fig. | for NML Cyg 




Fig. 13. Similar to Fig. | for U Her 



discussed in VO. The observed spectra are not directly pro- 
portional to /', as the minimum and maximum are closer 
together than predicted. We can express this narrowing as 
/n — AwL/Aw,nm, whcrc Aul is the FWHM of the maser 
feature and Aumm the separation between the minimum 
and maximum of the V-spectrum. Avmm is determined af- 
ter the spectrum is forced to be anti-symmetric by fitting 
the scaled down total power spectrum. In our spectra we 
find on average ~ 0.64, but sometimes as low as 0.45, 



which cannot be explained by LTE Zeeman models. In 
that case, when V oc /', we find /„ ~ 0.85. The non-LTE 
models partly explains this narrowing. As a result of the 
narrowing and rebroadening of the maser lines and of the 
hyperfine interaction, /„ depends on the maser brightness 
as shown in Fig. 
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The narrowing fraction decreases to 
w 0.6 in the case of highly saturated masers. However, 
as discussed above, the maser are indicated to be mostly 
unsaturated and thus the reason for the narrowing is still 
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Fig. 14. Fractional narrowing of the non-LTE models as 
a function of emerging brightness temperature. The solid 
line indicates the uni-directional maser, the dotted line is 
the bi-directional maser. Both models have wthormai = 1-0 
km/s. The long-dashed line is the theoretical value for the 
LTE Zeeman analysis, while the dashed-dotted line is our 
measured average narrowing value. 

not completely clear. We have performed some tests using 
bi-directional masers. As seen in Fig. |l^, these show more 
narrowing over the entire range, while the limit for the 
saturated masers is similar to the uni-directional masers. 
Thus, although the predicted narrowing fraction is still 
higher than observed, this is probably due to the assump- 
tion of linear masers in the models. Aside from the narrow- 
ing of the circular polarization profiles for bi-directional 
masers, the resulting A-p^-pi coefficients and correspond- 
ing magnetic fields do not change. We have used the uni- 
directional models for the fits, since the computation time 
increased considerably when calculating the bi-directional 
case. 

The above indicates that the narrowing of the V- 
spectrum is an important diagnostic for distinguishing 
which method yields the most reliable magnetic field 
strengths. A careful analysis of the total power intensity 
profiles will also be able to distinguish between different 
interpretations. The fits presented for our sources show 
that the non-LTE models do not fit the total intensity 
profiles as well as the fit of 3 hyperfine components for 
the LTE Zeeman models. However, this is due to the fact 
that it was impossible to completely explore the full pa- 
rameter space for the models. Fortunately, this has only 
a small effect on the magnetic field strengths inferred by 
this method, which has been included in the quoted errors. 
We thus conclude that the best estimates for the magnetic 



field strengths in the H2O maser region are those obtained 
using the non-LTE models. 

We have been unable to detect any linear polariza- 
tion for the sources in our sample. The upper limits are 
between 0.01 and 1.0%. This is consistent with previous 
observations by Barvainis & Deguchi (1989), and is a sep- 
arate indication that the masers are not strongly satu- 
rated. The absence of linear polarization excludes the non- 
Zeeman interpretation. We are therefore confident that 
the observed circular polarization is the signature of the 
actual Zeeman splitting and that our inferred values are 
the best estimates for the magnetic field. 

5.2. Magnetic fields in CSE 

We can compare our results with the magnetic fields mea- 
sured on the other maser species in the circumstellar en- 
velope. As remarked before there have been several obser- 
vations to determine the magnetic fields in the SiO maser 
regions (e.g. Barvainis et al. 1987, Kemball & Diamond 

1997) . Assuming that the circular polarization observed 
on these sources is caused by Zeeman splitting, fields be- 
tween 50 and 100 G have been inferred. This had already 
been predicted earlier from observations of the Zeeman 
splitting in OH masers, which indicated fields of 1 — 2 mG 
in Mira stars and up to « 10 mG in supergiants (Reid et 
al. 1979, Claussen & Fix 1982). Various other observations 
of OH masers have since then confirmed these milliGauss 
fields. The high field strengths on the SiO masers however, 
remain the topic of debate. It has been shown that fields 
of not more than w 30 mG can cause the high circular 
polarization observed; as described in the non-Zeeman ef- 
fect above. Also, Elitzur (1996) argues that the magnetic 
fields inferred on both the SiO and the OH maser could 
be lower by a factor of 8 - 10. 

Additionally, fields of tens to hundreds of Gauss in 
the SiO maser region would indicate fields of the order 
of 10^ G on the surface of the star. It was argued that 
such fields could not be produced by AGB stars (Soker & 
Harpaz, 1992). The high magnetic fields were determined 
by assuming the field strength varies with distance from 
the star by the relation B oc r". The exponent a depends 
on the structure of the magnetic field in the circumstellar 
envelope. A solar-type magnetic field has a = —2, while 
for a dipole medium a — —3. 

For a comparison between the different maser species, 
the results from this paper are shown in Fig. combined 
with results obtained on our sample of stars for the other 
masers. The figure displays the magnetic field strength as 
a function of distance from the star. We have indicated 
the actual measurements for the magnetic field with dif- 
ferent symbols for each star. The boxes indicate estimates 
for both magnetic fields and distance of the maser species 
based on observations of different stars (Johnston et al. 
1985; Lane et al. 1987; Barvainis et al. 1987; Lewis et al. 

1998) . We have separated the Mira variable stars and the 
supergiants. A major uncertainty is the distance of the 
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maser features for which the magnetic fields were deter- 
mined. The SiO masers are expected to exist in a nar- 
row region at approximately 1 stellar radius from the star 
(Reid & Menten 1997). The radial distances are obtained 
from the observation of their characteristic ring structure 
(e.g. Diamond ct al. 1994). Several observations have de- 
termined the extent of the OH maser regions of our sources 
(Diamond et al. 1987, Chapman & Cohen 1986). However, 
because the OH maser is strongly radially beamed these 
distances are only lower limits to the actual extent, which 
for the supergiant stars could be a few times 10^^ cm. 
The plotted measurements use the literature values quoted 
above as the value of the radial distance. The extent of 
the H2O maser region has also been the object of several 
studies (e.g. Bowers ct al. 1989, Lane et al. 1987). For 
our observations we have used the values quoted above. 
As the figure shows, our observations lie systematically 
above both the solar-type and the dipole magnetic field 
approximation. This can be caused by a selection effect. 
Our observations are heavily biased for the highest mag- 
netic fields, which occur closer to the star. Thus, while 
the distance values used in the plot are the outer edges of 
the H2O maser region, we most probably probe the inner 
regions. So using the solar-type magnetic field, we can es- 
timate the thickness of the H2O maser region. This ranges 
from 1 — 2-10^^ cm for the supergiant stars to « 2-10^* cm 
for U Her. This agrees with values found using the pump- 
ing models for H2O masers in late- type stars presented 
by Cooke & Elitzur (1985), and also with the estimates 
presented for S Per in Diamond et al.(1987). 

Another possibility is the existence of high density 
H2O maser clumps with frozen in magnetic field lines. 
For frozen magnetic fields the field strength varies with 
number density as _B oc n*^, with 1/3 < k < 1/2 from 
theoretical predictions (e.g. Mouschovias, 1987), where n 
is the number density. Density ratios between the maser 
clumps and their surrounding medium need to be between 
w 10 and « 1000 to explain the observations if we assume 
that the magnetic field strength in the medium surround- 
ing the clumps is similar to the field strength observed on 
the OH masers. While the higher values arc unlikely, the 
lower values are certainly reasonable. Because the obser- 
vations are most sensitive for the high magnetic fields, we 
are in this interpretation biased toward the highest dc;nsity 
maser clumps. A combination of the two effects described 
is likely. 

The pressure of the magnetic fields measured in the 
H2O maser region dominate the thermal pressure of the 
gas. The ratio between thermal and magnetic pressure is 
given by /3 = {SirnnkT) / . Here k is the Boltzman con- 
stant. Assuming a gas density of nu = 10^ cm""^ and a 
temperature of T « 10^ K at the inner edge of the H2O 
maser region, and taking B w 250 mG, /3 « 0.0-5, indi- 
cating that the magnetic pressure dominates the thermal 
pressure by a factor ?s 20. 

The inferred values for the magnetic field strength on 
the H2O masers are intermediate between the values de- 
termined on the OH masers and those determined for the 




log(R) [cm] 



Fig. 15. Magnetic Field strength, B, as function of dis- 
tance, R, from the center of the star. Dashed-dotted boxes 
are the estimates for the Mira variable stars, solid boxes 
are those for the supergiant stars. Similarly, the solid and 
the dashed-dotted line corresponds to a solar-type mag- 
netic field for the supergiants and the Mira variable stars 
respectively. The dotted line indicated a dipole magnetic 
field for supergiant stars. Symbols indicate observations. 
U Her is represented by triangles, S Per by the square, 
VY CMa by the crosses and NML Cyg by the hexagonal 
symbols. The dashed lines represent an estimate of the 
stellar radius. 

SiO masers with the standard Zeeman interpretation. This 
can be viewed as evidence that the measured circular po- 
larization on the SiO masers is caused by actual Zeeman 
splitting and that the non-Zeeman interpretation is un- 
likely even though the SiO masers exhibit up to 100% lin- 
ear polarization. Our observations agree with a power-law 
dependence of the magnetic field on distance to the star, 
of which the solar-type field seems most likely. An order of 
magnitude estimate for the magnetic fields on the surface 
of the star then gives fields of w 10^ G for Mira variable 
stars and « 10'^ for supergiants. Recently, Blackman et al. 
(2001) have shown that an AGB star can indeed produce 
such strong magnetic fields, originating from a dynamo at 
the interface between the rapidly rotating core and the 
more slowly rotating stellar envelope. 

6. Conclusions 

The analysis of the circular polarization of H2O masers de- 
pends strongly on the models used. The interpretation of 
our observations has shown that the non-LTE, full radia- 
tive transfer models presented by NW improve the accu- 
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racy of the fits. Some aspects of the observations however, 
remain unclear, even though the non-LTE approach can 
explain most of the narrowing of the circular polarization 
spectra. Multi-dimensional models seem promising to ex- 
plain the narrowing observed. We have found magnetic 
field strengths of « 150 — 500 mG for the supergiants and 
« 1.5 G for the Mira variable U Her. They occur in the 
most dense H2O maser spots at the inner boundaries of 
the mascr region and indicate strong magnetic fields at 
the surface of the star. The fields dominate the thermal 
pressure of the gas and are strong enough the drive the 
stellar outflows and shape the stellar winds. 
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